Cohesion-decohesion asymmetry in geckos 
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Lizards and insects can strongly attach to walls and then detach applying negligible additional forces. 
We propose a simple mechanical model of this phenomenon which implies active muscle control. We 
show that the detachment force may depend not only on the properties of the adhesive units, but also 
on the elastic interaction among these units. By regulating the scale of such cooperative interaction, 
the organism can actively switch between two modes of adhesion: delocalized (pull off) and localized 
(peeling). 

PACS numbers: 68.35.Np, 87.17.Rt, 87.85.gj, 87.10.Pq., 87.10. Hk, 62.20.Qp, 71.15.Nc 



I. INTRODUCTION 

Active mechanisms involved in biological adhesion in 
living systems are of broad theoretical interest in view of 
potential applications in bio-inspired adhesion devices. 
One of the most challenging issues concerns the reconcil- 
iation of strong adhesion pQ with easy detachment O [3J . 

Experimental studies reveal that biological adhesion 
at the organismic level is typically mediated by fibrillar 
microstructures which ensure a molecular level contact 
in the presence of surface roughness. The necessity of 
avoiding clustering entails a hierarchy of adhesion de- 
vices spanning a wide range of scales [4j [5] . In the case 
of geckos the macroscopic adhesion force at the level of a 
foot results from smaller forces at the scale of individual 
pads, each composed of tens of lamellae which in turn 
incorporate hundreds of thousands of setae. Each seta 
is split into hundreds of spatula shafts ending with spat- 
ula pads at the sub-micrometer scale where the adhesive 
force is ultimately generated by van der Waals forces [3] . 

Important insights into the functioning of microfibril- 
lar adhesive devices in geckos were obtained from AFM 
(Atomic Force Microscopy) experiments at the scale of 
spatulas [3] and setae [7J [5] and from attempts to ar- 
tificially microfabricate fibrillar microstructures [S]. In 
particular, these studies revealed a strong dependence 
of the adhesive forces on the angle formed by the se- 
tal shafts with the adhesion surface, which suggests that 
easy detachment may be achieved by active reorientation 
of the single microscopic fibers. Theoretical understand- 
ing of the angle dependence of adhesion at the level of a 
single spatula is mostly based on the study of Kendall's 
model 2J which has been recently generalized to account 
for the asymmetric attachment-detachment behavior of a 
single seta [TU] and for fiber tilting [TT] . The main idea is 
that decohesion can be described as peeling which implies 
that a Griffith's fracture takes place in an infinitely local- 
ized tip of a steadily propagating crack. An alternative 
model suggested in [7J [5] links the cohesion-decohesion 
asymmetry with a dependence of the cohesive strength 
on the tangential component of the force. An interesting 
attempt to reconcile such friction-based approach with 
Kendall's fracture model was proposed in |T5) where a fi- 



nite prestretch in the adhering layer was used to control 
the critical angle of adhesion. 

While fiber reorientation is clearly important for gecko 
adhesion, we propose in this paper a complementary 
mechanism that can be broadly characterized as the pos- 
sibility of active switching between localized (peeling) 
and delocalized (pulling off) fracture. We argue that 
the organism can control the modality of detachment by 
changing the level of coupling among individual fibrillar 
agents [13 . In contrast to Kendall's model, based on the 
assumption that the adhering layer can support only in- 
plane forces, we assume that this layer can have a shear 
stiffness which mimics bending resistance and is respon- 
sible for the cooperative effects. We also assume that the 
organism can actively switch between the regimes of high 
and low stiffness depending on the force that has to be 
exerted. We discuss potential mechanisms of how such 
control at different scales of the fibrillar micro-structure 
can be achieved and propose a strategy of gecko advance. 
Based on the experimental scaling relations, we conjec- 
ture that the same adhesion mechanism is operative at 
every scale of the hierarchy and propose a simple model 
justifying the observed power law force-length relations. 



II. DISCRETE MODEL 

To describe an elemental fibrillar adhesion layer we use 
the minimal Bishop-Peyrard (BP) model [THll6| . Con- 
sider a chain with n + 1 particles interacting through n 
linear springs and bound to a rigid substrate by break- 
able elastic elements (see Fig{T]). Assume that the par- 
ticles can move only in the vertical direction and denote 
by Ui the displacement of the particle i. The coupling 
between individual adhesive devices is controlled by the 
elastic energy 

Mb) = \g&1 

where G is the (shear) stiffness, Si — (u^+i — Ui)/l is the 
strain and I is the spring length. The cohesive energy at 
the micro scale can be modeled by a piece- wise quadratic 
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where k is the extensional stiffness, 7 is the adhesion 
energy density, b is the out of plane spatial scale and 
U r = v/276/fe is the limit displacement. 



n+1 



FIG. 1: Bishop-Peyrard model of decohesion. 

We assume that the organism applies to the pad a lo- 
calized loading. Denote by d the assigned displacement 
at the point n = 1 and assume that all other points are 
unloaded. To find the equilibrium force-displacement re- 
lation f{d) we minimize the energy 
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The corresponding energy landscape is complex |16j . 
however, we are interested only in a set of local mini- 
mizers that can be parametrized by the position of the 
decohesion front £ € [0, n], 
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Here the force is implicitly given by 
2nv 2 



, i = i + 1, n + l. 



(3) 
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(2f - 1 + coth % coth[(n + 1 - i)r\\) L 
where ~q is a solution of 

1 + l 2 /{2v 2 ) = coshM 

and 

V=y[Gjk 

is the internal length scale, defining the size of the cohe- 
sive zone. 




In Fig{2] we show the metastable branches f(d;£) 
parametrized by £; each branch ends at satisfy- 
ing = u r . If the dynamics is overdamped 
and the driving d(t) is quasistatic we obtain the loading- 
unloading hysteresis indicated in Fig{2ja) by bold lines. 
An alternative, hysteresis- free path of global energy mini- 
mization implied in Kendall's model is shown in Figj2jb) . 
Observe that in both cases the cohesion force exhibits a 
characteristic plateau with which we can associate the 
maximal force f m . The dependence of this threshold 
on the stiffness of the pad G plays the main role in 
the proposed mechanism. The quasistatic assumption, 



which this model shares with Kendall's model, is sup- 
ported by experimental observations that the attach- 
ment/detachment rates are independent of the geckos 
speed [IT]- 



III. CONTINUUM LIMIT 

A simple analytical expression for the function / m (G) 
can be obtained in the continuum limit, which may be a 
poor approximation at the scale of the whole foot (n ~ 
5), but turns out to be fully adequate at the level of 
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FIG. 3: Stiffness dependence of the adhesion threshold 
fm(G) for a system with u r = 1 mm, k = 1 MPa and 
L = 10 mm. Insets: crack configurations illustrating local- 
ized {L/v = 100) and diffuse (Ljv — 1) cohesive zones. 



the setae (n ~ 50). To study this limit we fix the total 
length L = nl and assume that n — > oo and I —> 0. The 
continuum energy takes the form 



$(u) 



(5) 



and its minimization is straightforward |161 118] . Given 
the loading u(0) = d we obtain 



u ( (x)- 
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where C is the detached fraction of the pad which 
parametrizes the relation between the boundary displace- 
ment 

d = u r + (L(f/G) 

and the applied force 

/ = v /2 7 6Gtanh[ J L(l - Q/v\. 

The maximum value of the force f m , representing the 
detachment threshold, is attained at ( = 1 and d = u r : 



f m {G) = v / 2jbGtanh(LVk/VG). 



(7) 



In the limit when the external length scale is much larger 
than the internal length scale (L ^> v) we get the asymp- 
totics 



f m ~ 

In the opposite case we obtain 

f m ~ ^bkL. 

The value of the critical force and the structure of the 
associated cohesive layers is shown in Fig(3] Notice that 
for L»v the crack has a narrow tip and this regime can 
be associated with the peeling mode of fracture. Instead 



for L <C v the cohesive zone spreads along the whole pad 
and this regime can be associated with a pull off mode. 
Interestingly, for L S> v the area under the hysteresis 
loop 

Q = 2-ybL 

representing the detachment energy is exactly twice as 
big as in the case when L <C v, high dissipation at small 
force is due to much larger displacement. In this sense 
transition from pull off to peeling is similar to the tran- 
sition from brittle to ductile fracture. 

Our main assumption in what follows is that the stiff- 
ness of the linear springs, mimicking the stiffness of the 
gecko's pad, can be actively varied by the organism. The 
feasibility of such control is clear from the fact that the 
gecko rolls in for attachment (active shortening and thick- 
ening of the digits) and rolls out for detachment (active 
lengthening and thinning of the digits). It is also known 
that digital hyperextension anticipates each attaching 
and detaching event and that the musculo-tendinous sys- 
tem may influence single lamellae in the process of con- 
trolled rolling [T31 HZ]. Experiments with geckos [T5] . 
frogs [IH] , and ants |2D] also show that at the macroscale 
the easy release is achieved through the localization of 
the cohesive region. Yet another argument in support 
of the active muscular control comes from the observa- 
tion [5T] that to simplify horizontal walking geckos keep 
the hyperextended state and start activating the attach- 
ment mechanism only at sufficiently high slope requiring 
stronger adhesion. 

To get a rough estimate of the required stiffness vari- 
ation we observe that the adhesion force in geckoes is 
about one order of magnitude less than the detachment 
threshold [35]. According to ^ this corresponds to two 
orders of magnitude in stiffness variation, which is com- 
patible with the data on geckos forced to detach [23] . 
Different physical mechanisms may be employed to regu- 
late the coupling among adhesive fibrils at different struc- 
tural levels. Thus, at the cellular level the dynamic fil- 
amentous actin network is known to be very soft at low 
stress, but can stiffen up to three orders of magnitude in 
response to stresses [2"lH2"o] . Such stresses can be gen- 
erated internally through molecular motors; moreover, 
constant remodeling allows the cytoskeleton to remain 
in a marginally stable state and easily switch between 
softening and stiffening [3Z] HE] • At elevated stresses a 
quick transition to softening may also be related to the 
unfolding of the cross-linkers such as filamin [29] . Within 
muscle sarcomeres the effective stiffness can vary with the 
number of myosins attached to actin fibers and will also 
be affected by the unfolding-refolding transition in titin 
[30] , Notice also that parameter v controlling the mode 
of detachment can vary not only because of the stiffness 
G but also because of the stiffness k. The latter depends 
on the aspect ratio of the adhesive elements and may be 
controlled by the organism through capillarity induced 
self-assembly, modulated by secretion or evaporation of 
liquids responsible for the interaction between the fibrils 
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At the macroscale, the variability of shear modulus 
may be associated with the reversible development of mi- 
cro defects inside the tissue architecture with subsequent 
internal healing of this damage through remodeling. Fol- 
lowing a classical approach in damage mechanics [33|,l34"]. 
we can introduce an internal variable a, with a = repre- 
senting the undamaged state (stiffest configuration) and 
a = 1 - the damage-saturated state (most compliant con- 
figuration). Then G = G(a) with G'(a) < 0. Since the 
detachment force depends on the level of damage 

fm(a) = fm(G(a)), 

while the critical displacement does not (d m = u r ), we 
can write the effective elastic stiffness for the attached 
state in the form 



E(a) = \jkG{a)tanh[LVk/G{a)}. 
Observe that E'(a) < 0. 



/(N) 




FIG. 4: Schematic representation of the proposed 

attachment-detachment strategy (see text). The inset shows a 
generic force-displacement path associated with stiffness vari- 
ation AB and the decomposition of the corresponding external 
work W into the elastic energy <3> and the dissipated energy 
Q. 



IV. ATTACHMNT-DETACHMENT CYCLE 

We can now propose a mechanical strategy of geckos 
advance. Assume that only one scale of the fibrillar mi- 
crostructure is involved, that the response is quasistatic 
and overdamped, and that the continuum limit is valid. 
To fix ideas, consider ceiling walking and suppose that 
the cycle of attachment-detachment starts at an attached 
state A (see inset in Fig. [4| where the pad is in the 
stiff configuration G(0) associated with the high thresh- 
old jm(0). To attain the detachment below this thresh- 
old, gecko can induce a reversible "damage process" de- 
scribed by a force-displacement relation / = f(d) where 
d = d(a) and E(a) = f (d(a)) / d(a) . In the process of 



stiffness variation the external work 

pd(u) 

W = / f(d) dd>0 

JP/E(0) 

(done by gravity) is partially stored as elastic energy 

l<ly\=P 2 /(2E(0))-E 2 (a)d(a)/2 

and the rest is dissipated into heat 
f«l, 

Q = - ~E'(a)d 2 (a)da. 
Jo 2 

To minimize dissipation while maintaining a stable me- 
chanical response with f'(d) > 0, the animal must ensure 
that damage advances at constant force f(d) = P (AB in 
Figj4]). In this case exactly half of the work is dissipated 
and we can write 

\$]=Q = (P/2)(u r -P/E(0))>0. 

Damage at constant stress has been observed in many 
rubber-like materials and linked to inherent energy non- 
convexity [35] : reversible structural changes at constant 
stress are also characteristic for muscle tetanus [55] . 

After the critical force is decreased, the pad can be 
pulled away by peeling (path BC in Figj4]) . The detach- 
ment ends with an abrupt decohesion at point C in Figj4] 
In order to reattach, the gecko can follow a reverse path 
EF. As the foot is placed on the surface, the displacement 
d is gradually decreased and the attachment takes place 
at the point F through "inverse peeling" at d = u r . An 
instantaneous force jump brings the system to the point 
G and allows the animal to place some weight on the 
foot. To secure a robust attachment, the gecko must re- 
verse the damage and induce active stiffening (trajectory 
GH). 

The system heals the damage by remodeling the dam- 
aged configuration with a = 1 back into the virgin con- 
figuration with a — while increasing stiffness and de- 
creasing displacement. This requires work which is now 
done by the gecko. The energy comes from metabolic 
sources M < and is released due to elastic unloading 

[$] = (/ m (l)/2)(/ m (l)/£(0) - u r ) < 0. 

If we again assume that healing (remodeling) takes place 
at constant stress (minimum metabolic energy path), 
we obtain [<i>] = M. Experiments show that geckos in 
the compliant state have a very low detaching thresh- 
old /m(l) [37] which means that the metabolic energy 
required for such stiffness increase is also low. 

After the state of high stiffness is reached (point H in 
Fig(4]) the peeling mode is deactivated because the force 
required for detachment / m (0) is now large. Therefore 
more weight can be shifted to this foot (path HI) and 
another foot can undergo the detaching-attaching cycle. 
Overall, the detachment process BCD with decreasing 
force on the two detaching feet must take place simul- 
taneously with the attachment process FGHA involving 
the other two reattaching feet [17] , 
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FIG. 5: Scale free nature of the adhesion mechanism. Solid 
line interpolates the experimental values of the adhesion forces 
[H H IHEl [SB [23] E3 |39]; dashed line, of the friction forces 
|37| . Insets show nested computations based on the BP model 
with adhesion links behaving as shown in the right upper 
corner. At the smallest scale of the elemental fibrilla u r = 200 
nm, L — 10 3 nm and k = lOOMPa. The parameters u e , 
u r , and k at larger scales have been computed iteratively. 
Other parameters: fibrilla — > seta, n = 50, L = 100 nm, 
G = 0.2 mN; seta -> lamella, n = 50, L = 0.5 mm, G = 1 N, 
lamella — > pad, n = 20, L — 3 mm, G = 300 N, pad — > foot: 
n = 5, L = 20 mm, G = 20 N. 

V. HIERARCHICAL ARCHITECTURE 

To understand the role played by the hierarchical 
micro-fibrillar architecture of the adhering pad [T7J1 [3H] 
we compare adhesion forces at different spatial scales 
(see Fig|5|. Using the experimental data [TJ [3] EJ- 
[Hl EH I23[ 1371 155] we may deduce a scaling relation 

f ~lP 

with exponent /3 = 2 ± 0.5. In Figj5]we show the power 
law together with the scaling of friction forces reported 
in |37| ; our results are consistent with the observation 
that the friction force is usually five times higher than 
the adhesion force |2"3] . 

While the exact origin of these empirical power law re- 
lations is not known, we can propose the following plausi- 
ble explanation. Following [40] we suppose that the archi- 
tecture of the fibrillar adhesive system is designed to en- 
sure that in contact with rough surfaces the fibrils buckle 
simultaneously at all scales to ensure maximum folding 



which accommodates fractal roughness. If at hierarchi- 
cal level i, with spatial scale U, the system is marginally 
stable against buckling, then fa = cEl?, where E is the 
elastic modulus and the constant c depends on the shape 
of the cross section, aspect ratio, and the boundary con- 
straints. We assume the simplest allometric law when 
both c and E are scale invariant. Now, if the fibrils at 
a finer level i cover the tips of the fibrils at the coarser 
level i + 1 densely, which is known as Leonardo's rule 
|42) . then = riilf where rii is the number of fibrils 
at the level i (see the scheme at the bottom right corner 
of Fig. [5|. Since = n t fi, we obtain = cEl? +1 , 
which means that the force is critical also at the i+1 level. 
In view of our neglect of collective modes of instability 
[43] . we can only tentatively conclude from this reason- 
ing that the observed scaling supports the idea that the 
whole structure can be marginalized simultaneously. An 
important feature of this scaling, however, is that stress 
is uniform which has been proposed previously as a cri- 
terion of optimality in several biological and engineering 
systems }44j . 

The power law scaling is indicative of a scale-free de- 
tachment mechanism. We can model it in our framework 
by the appropriate renormalization of the parameters in 
a scale-generic BP model shown in the upper corner of 
Fig(5] where u e is the elastic threshold and u r is the 
detachment displacement. The parameters u e , u r , and 
k can be computed at each scale iteratively by using a 
series of nested BP models, whereas the parameter G 
characterizing the elastic coupling can be chosen at each 
scale to match the experimentally measured detachment 
threshold. 

At the smallest scale of a spatula the adhesive proper- 
ties of the fibrils can be described by the energy density 
([l]) with parameters u r and k available from experiment 
[Sj. The behavior at the next scale (setae) can be ob- 
tained numerically by simulating an overdamped gradi- 
ent flow dynamics for a quasistatically driven BP system 
with n — 50 spatulae; the value of G is chosen to ensure 
the maximum measured adhesion force of 40/iN [7] . The 
overall behavior at this scale matches the experimental 
results in [7] showing an elastoplastic range ending with 
an abrupt detachment. At the level of a lamella, we con- 
sider n = 50 elastic-plastic elements with constitutive 
parameters obtained from the spatulae scale simulations. 
At the next level of a toe, we need to model lamellae 
with n — 20 and finally at the level of a foot n — 5 (toes) 
and the problem becomes strongly discrete. The results 
at this last level match observations showing digitized 
detachment of the toes [T3J • 



VI. CONCLUSIONS 

By studying a prototypical system, we have shown that 
both the force threshold and the dissipation associated 
with reversible adhesion can be modified by active con- 
trol of the coupling among individual adhering elements. 
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The possibility of the ensuing multi-path adhesion [IS] PRIN 2010-11: "Dinamica, stabilita e controllo di strut- 
reconciles strong binding with easy debinding, which is ture flessibili" . 
at the base of the observed agility of lizards and insects 
running on inclined surfaces. The proposed mechanism 
has a scale- free hierarchical structure, which is typical for 
biological systems at all levels of organization. 
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